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Quantum  Dynamical  Model  of  Laser/Surface- Induced 
Predissociation:  Multi photon-Multiphonon 
-  Processes  and  Photon/Phonon-Dressed  States* 

Jui-teng  Lin  and  Thomas  F.  George + 

Department  of  Chemistry,  University  of  Rochester 
Rochester,  New  York  14627  USA 

By  the  dynamical  energy  populations  of  a  photon/phonon-dressed 
quantum  system,  laser-stimulated  surface  predissociation  (LSSPD)  of 
an  adsorbed  diatomic  molecule,  characterized  by  the  laser  pumping 
rates,  the  phonon-induced  damping  factor  (level  width),  the  field 
detuning,  and  the  phase  detuning  (Landau- Zener  form) ,  is  studied. 

The  optimal  detunings  for  the  enhancement  of  LSSPD  are  shown  by  the 
numerical  results  of  a  set  of  coupled  generalized  quantum  equations 
(GQE) .  The  important  features  of  LSSPD  and  the  possible  predissocia¬ 
tion  of  a  multilevel  multistate  system  via  multiphoton-multiphonon 
processes  are  discussed. 
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1.  Introduction 

The  interaction  of  laser  radiation  with  matter  (atoms  and 
molecules)  of  homogeneous  systems/  e.g. ,  multiphoton  dissociation 
and  ionization,  have  been  intensely  and  systematically  studied  since 
the  tunable  and  high-power  lasers  arrived  at  the  beginning  of  the 
1970's.  The  study  of  the  influence  of  lasers  on  the  dynamics  of 
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heterogeneous  gas/surface  systems,  however,  is  still  in  its  infancy. 

The  complexity  of  laser-stimulated  surface  processes  (LSSP)  is  not  only 
due  to  the  heterogeneity  of  the  system  (where  the  many-body  effects 
of  the  surface  atoms  are  essential)  but  also  due  to  the  laser-stimulated 
(or  enhanced)  dynamical  phenomena  of  the  adspecies,  e.g.,  surface 
migration,  desorption,  decomposition,  predissociation,  recombination 
and  other  surface  rate  processes. 

The  influence  of  IR  radiation  (with  vibrational  excitations) 

in  LSSP  has  been  theoretically  studied  in  our  laboratory.  Some  of 

the  important  features  and  the  related  methodology  in  our  previous 
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model  systems  are: 

(1)  The  average  excitation  (photon  energy  absorbed)  of  the  ac¬ 
tive  mode,  the  absorption  cross  section  (or  the  line  shape)  and  the 
pumping  rate  of  the  radiation  are  characterized  by  the  field  intensity 
and  detuning,  the  T^  and  T2  relaxation  factors,  the  phonon  coupling 
factor  and  nonlinear  effects  due  to  the  anharmonicity .  The  results 

have  been  obtained  either  by  a  quantum  approach  (in  the  Heisenberg- 

6“  8 

Markovian  picture)  or  by  a  stochastic  classical  approach  via  the 
generalized  Langevin  theory. ^  ^ 


(2)  The  mechanism  of  laser-controlled  surface  phenomena  (selec¬ 
tive  desorption  and  migration)  may  be  investigated  by  combining  the 
Langmuir  kinetic  equation  with  the  laser  rate  equation  (in  the  energy 
space)  or  generalized  master  equation  (in  the  energy  plus  lattice 
site  space) . 5 

(3)  By  using  Heisenberg  equations  of  motion  (in  the  rotating- 
wave  approximation) ,  we  are  able  to  find  the  dynamical  level  popula¬ 
tions  of  different  normal  modes  and  investigate  the  nature  of  LSSP, 
mode-selective,  local  heating  (molecular-selective)  or  simply  thermal 

q  in 

(nonselective)  heating  of  the  system.  ' 

In  this  paper,  we  shall  study  another  aspect  of  LSSP,  namely 
laser-stimulated  surface  predissociation  (LSSPD) ,  by  investigating 
the  synergistic  effects  of  surface  and  two  lasers,  providing  both 
vibrational  and  electronic  transitions  of  the  adspecies,  which  may 

dissociate  only  with  negligible  probability  (at  low  temperature)  in 
the  presence  of  the  surface  or  a  single  laser.  Since  the  phonon- 
induced  damping  factor  and  the  surface  potentials  of  the  adspecies 
are  not  experimentally  available  to  date,  instead  of  focusing  on  a 
specific  real  system  we  shall  study  the  general  dynamical  feature 
of  LSSPD  by  using  parameterized  quantities.  The  main  purpose  of  the 
present  paper,  therefore,  is  to  predict  and  suggest  some  possible 
enhancement  mechanism  of  LSSPD  for  real  systems  based  on  the  result¬ 
ing  feature  of  a  model  system. 
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2.  LSSPD  of  Photon/Phonon-Dressed  Multilevel  System 

We  consider  a  heterogeneous  model  system  consisting  of  an  optically 
active  mode  and  a  group  of  inactive  modes  (phonon  modes  plus  the  in¬ 
active  modes  of  the  adspecies)  and  subject  to  two  lasers  of  different 
frequencies.  While  the  active  mode  is  vibrationally  excited  by  the 
first  laser  (IR) ,  the  adspecies  is  also  electronically  excited  by 
the  second  laser  (UV  or  visible) ,  where  the  curve  crossings  among  the 
electronic  states  in  turn  cause  the  vibrational  predissociation  of 
the  adspecies  (diatomic  molecule)  on  the  solid  surface.  When  the 
species  is  adsorbed  (particularly  in  the  case  oi  chemisorption)  on 
a  solid  surface,  its  vibrational  levels  may  be  shifted  and  broadened 
not  only  by  its  intramolecular  anharmonic  coupling  but  also  by  the 

phonon  coupling  (usually,  multiphonon  processes)  and  the  dynamical 

£ 

coupling,  e.g.,  surface  migration-induced  or  collision-induced. 

The  time-dependent  Schrddinger  equation  of  the  photon/phonon- 
dressed  state  may  be  written  as 


Hq  is  the  effective  field-free  Hamiltonian  of  the  phonon-dressed 
state  characterized  by  a  complex  form 


*C-H- 


where  H_  is  the  unperturbed  molecular  Hamiltonian,  and  <Su)  and  y 

u  mn  mn 

are,  respectively,  the  frequency-shift  and  the  level-broadening 
(damping  factor)  of  the  phonon-dressed  state  |m,n>  [m  and  n  represent 
the  m-th  electronic  and  the  n-th  vibrational  states,  respectively]. 
Considering  only  the  dominant  factor  of  the  frequency-shift  and  the 
level-broadening,  i.e.,  via  multiphonon  coupling,  these  two  phonon 
parameters  may  be  expressed  by®'^'^ 

IKcpfit  , 

-n  —  (iT/ka^)  [  -Mf  (  1  ] 

where  K  is  the  coupling  factor;  p  and  are  the  density  of  states 

and  frequencies  of  the  phonon  modes,  respectively;  p  is  the  smallest 

order  of  the  multiphonon  coupling  given  by  the  conservation  of  energy 

p  *  u>  /u>  ;  a)  and  co  are  the  Debye  frequency  and  the  frequency  of 
mn  D  D  mn  J 

the  active  mode  (for  a  transition  between  states  |m,n+l>  and  |m,n>) , 
respectively;  and  T  is  the  surface  temperature. 

The  interaction  Hamiltonian  which  couples  the  active  mode  with 
the  lasers  is  given  by 

H'(t )  =  V,  +  Vz  (4) 

where  «  u^E1  (due  to  the  first  laser)  is  proportional  to  the  deriva¬ 
tive  of  the  dipole  moment,  u£,  (with  respect  to  the  active-mode  normal 


(3a) 

(3b) 

(3c) 


k 
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coordinate)  and  =  uE2  {due  to  the  second  laser)  is  proportional 
to  the  dipole  moment  of  the  excited  species,  u.  and  E2  are  the 

envelopes  of  the  electric  field  (assumed  to  be  slowly  varying  or  con¬ 
stant  in  time)  of  the  first  laser  for  vibrational  excitation  (with 
frequency  w^)  and  of  the  second  laser  for  electronic  excitation  (with 
frequency  w2)  ,  respectively. 

Using  this  photon  phonon-dressed  state  governed  by  a  complex 
effective  Hamiltonian,  we  may  then  introduce  the  generalized  quantum 
equations  of  motion  (GQE)  in  the  Heisenberg-Markovian  picture  and 
the  rotating-wave  approximation  (RWA) : 


where  W  (R(t) ]  is  the  real  part  of  the  diagonal  element  of  the  field- 
mn 

p  f  f 

free  phonon-dressed  effective  Hamiltonian  H^  is  the  basis-state 
representation  of  the  unperturbed  vibronic  state  4>^(r,R);  r  and  R 
represent  the  coordinates  of  the  electrons  and  of  the  nucleus,  re¬ 
spectively;  6  .  and  6  .  are  Kronecker  delta  functions  assigned  for 
r  J  m  m  n  m 

the  first  and  second  lasers  which  cause  vibrational  (m’=m)  and  elec¬ 
tronic  (n'=n)  transitions,  respectively. 

We  note  that  GQE  (eq.(5)]  reduces  to  the  usual  Heisenberg  equa¬ 
tions  of  motion  (in  RWA) ^  in  the  atomic  viewpoint,  i.e.,  where  one 


(S) 
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atom  is  viewed  as  perturbed  by  another  atom  and  W  (R)  is  a  complex 
time-independent  quantity,  W^fR)  =  Emn  ~  ^^mn^^ '  an<^  ky  using  the 
transformation  ( t)  =  (t)  exp  [y^t/2] , 

y*  ,ti 

Furthermore,  GQE  may  be  reduced  to  the  semiclassical  equations  of 
motion  (in  the  molecular  viewpoint)  for  phonon-free  case,^  i.e., 
Ymn=0 .  It  is  also  worth  noting  that  the  equations  of  motion  in  the 
molecular  viewpoint  may  be  developed  by  simply  including  the  secular 
terms  on  the  right-hand  side  of  eq. (7)  and  using  the  transformation 

Q «> ■  f fcj*’ dt‘\  . 


3.  LSSPD  of  a  Two-Level  Two-State  System 

GQE  shown  in  eq.(5),  describing  a  multilevel  (vibrational)  multi¬ 
state  (crossing  electronic  curves)  photon/phonon-dressed  system,  is 
in  general  characterized  by  multiphoton-multiphonon  processes.  These 

coupled  first-order  differential  equations  in  principle  may  be  de- 

1 8 

coupled  into  a  single  high-order  differential  equation,  or  solved 
by  the  Laplace  transform  technique .  ^  However,  the  analytic  solutions 
are  not  usually  available  due  to  the  vibrational  couple  terms  but 
also  due  to  the  time-dependent  classical  trajectory  factor  W^fRft)]. 

To  investigate  the  synergistic  effects  of  the  lasers  and  the  surface 
phonons  by  some  tractable  means,  we  shall  consider  the  physical  situa¬ 
tion  that  the  adspecies  may  be  predissociated  by  absorbing  two  photons 
(one  IR  and  one  UV,  or  visible  photon) .  We  shall  also  assume  that 


there  is  only  one  electronic  curve  crossing  (or  one  of  the  crossings 
contributes  the  dominate  predissociation  channels) . 


For  a  two-level  two-state  system,  GQE  [eq. (5) ]  reduces  to  a  set 
of  coupled  equations  of  the  amplitude  vectors  C^(t),  C^2(t)  and 
C2c(t)  representing  the  states  of  bound-ground,  bound-excited  and 
the  continuum  (predissociated)  states,  respectively: 

4 

Cn  (t >  ~  C]2  ,  (9a) 

c.  ft)  =  -a  v (*&)  C  ^  ’ 

12  r  \  //  (9b) 

C  it)  =  -i  G[  (i  cCt%  )  CfZ  . 

(9c) 


The  pumping  rates  v  and  G  are  given  by  [from  eq.(6.a)],  V=V^  « E^ 

2c 

and  G  =  V^2  aE^;  y  is  the  level  width  of  the  phonon- dressed  upper 
level  referred  to  as  the  damping  factor  of  the  pumped  system;  A  =  -u> 

\1S/ 

is  the  field  detuning  with  the  fundamental  frequency  of  the  phonon 

dressed  state;  and  a=  (S2-S^)Vg/1i  is  the  phase  detuning  obtained  by 

the  usual  Landau-Zener  model,  i.e.,  W2c  ( R) -W12  (R) -hco2*'at .  "*■ 7  There 

are  some  important  features  of  the  above  coupled  equations: 

(i)  The  two-photon  predissociation  processes  are  two-step  excitations 

governed  by  the  pumping  rates  V  (due  to  the  first  laser)  and  G  (due 

to  the  second  laser),  (ii)  due  to  the  damping  factor  y  (energy  sink), 

both  the  populations  of  the  bound  state  |c^(t)  |  and  |c^2(t)  |  are 

2 

decaying,  while  the  probability  of  predissociation  PD(t)  =  |C2c(t)  ! 
will  rise  to  a  saturated  value,  and  (iii)  the  oscillatory  behavior 
due  to  the  phase  detuning  [exp(iat  ) ]  is  significantly  different  from 


that  of  the  field  detuning  (exp(-iAt)],  and  we  shall  see  later  that 
the  corresponding  optimal  detunings  show  different  character. 


The  analytical  solutions  of  the  coupled  equations  are  available 
in  some  specific  situations  where  one  may  use  perturbation  methods. 
However,  we  are  more  interested  in  the  exact  solutions  which  are 
available  by  using  numerical  method,  for  e.g.,  Runge  Kutta.  The 
numerical  results  of  the  level  population  P^=  |c^2(t)|  (for  photon 
energy  deposited  in  the  active  mode) ,  the  predissociation  probability 
PD  =  1  an<3  t*ie  energy  population  of  the  phonon  (bath)  modes 

PB  =  1  -  |Cj_-j_(t)|  -  P^  -  PD  are  shown  in  Figs.  1-4.  There  are  five 

parameters  characterizing  the  dynamical  feature  of  LSSPD:  (V,G,y,A,a) 
For  systematic  analysis  of  the  numerical  results,  we  investigate  the 
following  physical  features: 

(i)  Fig.  1  shows  the  effect  of  the  phase  detuning  a  (for  y=0)  on  the 

steady-state  predissociation  probability  (Pp‘S‘).  It  is  seen  that 

instead  of  a=0,  the  optimal  phase  detunings  for  p^‘s*  having  peak 

values  are  a* =±1.0  and  a* =±3.5;  and  for  |a|>3.5,  the  values  of 
S  •  s 

PD*  '  show  exponential  decaying.  This  is  a  significantly  different 

feature  from  that  of  the  field  detuning,  where  there  are  no  dips  in 

the  absorption  profile  and  the  optimal  detunings  are  A*=0  for  a  two- 

level  system  while  A*>0  (red-shift)  for  a  multilevel  system  with 

6  7 

anharmonic  nonlinear  effects.  ' 

(ii)  The  effect  of  the  phonon-induced  damping  factor  (y)  is  shown  in 
Fig.  2 ( A) - (D) ,  for  (V,G,A,y , a)  =  ( 4 , 1 , A* , y , a*)  with  the  optimal 
detunings  A*=0  and  a*=3.5  and  the  damping  factor  increasing  from  y=0 


! 

i 


(surface  free)  to  y =  1  (strong  phonon  coupling). 

(iii)  The  effect  of  the  phase  detuning  on  the  dynamical  populations 
and  the  time-dependent  predissociation  is  shown  in  Fig. 3 (A) - (B) ,  for 
the  cases  of  exact  phase  resonance  (a=0)  and  far-off  resonance  (a=9) . 
We  note  that  P®,S"  =  0.06  for  a=0 ,  and  P®'s‘=  0.15  for  a=9 . 

(iv)  Finally,  we  show  the  high  Rabi  oscillation  cases  (governed  by 
the  first  laser  pumping  rate  V)  for  y=0  and  y=0.1  with  the  other  op¬ 
timal  phase  detuning  a*=l.  The  other  important  physical  quantity  is 

the  mean  predissociation  rate  which  may  be  approximately  defined  as 
—  s  s 

R  ~PD"  ‘/t*,  with  t*  being  the  rising  time  of  the  predissociation 
probability  P  (t).  We  can  estimate,  from  Fig. 2(C)  and  2(D),  that  the 
mean  predissociation  rates  (per  unit  time  of  1/G)  decrease  from 
R  -0.21  to  R  -0.14  when  the  damping  factor  increases  from  y=0 . 1  to 
7=1. 

4.  Discussion  and  Conclusion 


As  an  extension  of  the  model  system  described  in  the  previous 
section,  we  shall  now  suggest  a  possible  enhancement  of  surface  pre¬ 
dissociation  of  a  multilevel  (vibrational)  and  multistate  (electronic 
curve-crossing)  system  via  the  combination  of  two  lasers  as  follows: 

(1)  LSSPD  may  be  enhanced  by  absorbing  two  laser  photons  where 
the  optimal  detuning  condition  (A=A*  and  a=a*)  is  available  when  both 
lasers  (with  frequencies  and  w^)  are  tunable. 

(2)  Effective  predissociation  channels  induced  by  surface  mag¬ 
netic  field  can  be  opened  in  the  absorption  of  one  or  few  IR  photons, 

i.e.,  preparing  the  adspecies  in  a  vibrational  level  above  the  highest 

19 

surface-split  electronic  states  and  providing  the  external  photon 
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energy  source  to  balance  the  energy  loss  due  to  the  phonon  coupling 
of  the  adspecies.  We  note  that  in  this  case  the  net  enhancement  of 
LSSPD  is  governed  by  the  competition  between  the  surface-split  (which 
opens  more  possible  predissociation  channels)  and  the  surface-induced 
energy  damping.  Therefore,  the  gaining  of  a  significant  steady-state 

g  g  _ 

predissociation  PD‘  *  or  the  rate  of  predissociation  R  is  governed  by 
the  competition  between  the  (IR)  multiphoton  and  the  (surface)  multi- 

g 

phonon  processes,  and  moreover  by  the  coupling  strength  of  the  Landau- 

Zener  local  transition  probability  characterized  by  the  pumping  rate 

17 

of  the  second  laser  (G)  and  the  phase  detuning  (a) . 

(3)  For  a  multilevel-multistate  system  (see  Fig. 5),  the  enhance¬ 
ment  of  LSSPD  may  also  be  generated  by  multiphoton  processes  where 

the  average  excitation  of  the  multilevel  active  mode  (with  energy 

6  20 

population  governed  approximately  by  a  Poisson  distribution)  ' 

matches  one  of  the  peak  values  of  the  Franck-Condon  factors  which  show 

21 

structure  patterns  as  a  function  of  the  vibrational  level. 

In  conclusion,  we  have  shown,  by  the  numerical  solutions  of 
a  set  of  coupled  GQE,  how  the  synergistic  effect  of  two  lasers  and 
the  surface  phonons  on  the  dynamical  feature  of  LSSPD  are  character¬ 
ized  by  the  pumping  rates  (V,G) ,  the  multiphonon  coupling  factor 
(y) ,  the  field  detuning  (A)  and  the  phase  detuning  (a) .  For  sig¬ 
nificant  enhancement  of  LSSPD,  we  expect  a  small  ratio  of  y*v  and  match¬ 
ing  detunings  to  the  optimal  condition,  i.e.,  A=A*  and  a=a*. 

We  are  grateful  to  Drs .  M.  Hutchinson  and  K.  S.  Lam  for  their  helpful 
discussions . 
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Figure  Captions 

Fig.  1.  The  steady-state  predissociation  probability  (P®*s*)  as 
a  function  of  the  phase  detuning  (a)  for  (V,G,y,A)  =  (4, 1,0,0). 

V, G, y  and  A  are  the  pumping  rates,  phonon-induced  damping  factor 
and  the  field  detuning,  respectively. 

Fig.  2.  The  time  evolutions  of  the  energy  populations  of  the  active 
mode  (P  )  and  the  phonon (bath)  modes  (P  ) ,  and  the  predissociation 
probability  (PD)  for  (V,G,y,A,oi)  =  (A)  (4, 1,0, 0,0)  [phonon-free  at  exact 

resonance],  (B)  (4 , 1 , 0 , 0 , 3 . 5)  [phonon-free  at  optimal  detuning], 

(C)  (4, 1,0. 1,0, 3. 5)  [weak  phonon-coupling],  and  (D)  (4 , 1 , 1 , 0 , 3 . 5) 

[strong  phonon-coupling] . 

Fig.  3.  The  dynamical  feature  of  PA,  P0  and  PD  for  (V,G,Y,A,a)  = 

(A)  (4, 1,0. 5, 0,0)  [exact  resonance  with  phonon-coupled],  03(4,1,0.5,0,9) 
[far-off  phase  detuning] . 

Fig.  4.  The  dynamical  feature  of  P^,  Pg  and  PQ  for  the  cases  of  high 
Rabi  frequency  (pumping  rate  V)  with  (V,G,Y,A,a)=  (A)  (10,1,0.1,0,1), 

for  phonon-coupled  and  (B)  (10,1,0,0,1),  for  phonon-free  system. 

Fig.  5  (A)  Schematic  energy  diagrams  of  a  multilevel-multistate  system 
subject  to  two  lasers  with  frequencies  u>  and  w',  respectively. 

(B)  The  vibrational  energy  population  of  the  photon/phonon-dressed 
state  PR  vs.  n  (vibrational  quantum  number). 

(C)  The  structure  pattern  of  the  Franck-Condon  factors  for  a  multistate 
(curve-crossing)  system. 
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